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Graphene hosting a pair of collinear adatoms in the phantom atom configuration has density of 
states vanishing in the vicinity of the Dirac point which can be described in terms of the pseudogap 
scaling as cube of the energy, A oc |e| 3 which leads to the appearance of spin-degenerate bound 
states in the continuum (BICs) [Phys. Rev. B 92 , 045409 (2015)]. In the case when adatoms are 
locally coupled to a single carbon atom the pseudogap scales linearly with energy, which prevents 
the formation of BICs. Here, we explore the effects of non-local coupling characterized by the Fano 
factor of interference qo, tunable by changing the slope of the Dirac cones in the graphene band- 
structure. We demonstrate that three distinct regimes can be identified: i) for qo < q c i (critical 
point) a mixed pseudogap A oc |e|, |e| 2 appears yielding a phase with spin-degenerate BICs; ii) near 
q 0 = q c i when A oc |e| 2 the system undergoes a quantum phase transition (QPT) in which the 
new phase is characterized by magnetic BICs and iii) at a second critical value qo > q c 2 the cubic 
scaling of the pseudogap with energy A oc |e| 3 characteristic to the phantom atom configuration is 
restored and the phase with non-magnetic BICs is recovered. The phase with magnetic BICs can 
be described in terms of an effective intrinsic exchange field of ferromagnetic nature between the 
adatoms mediated by graphene monolayer. We thus propose a new type of QPT resulting from the 
competition between two ground states, respectively characterized by spin-degenerate and magnetic 
BICs. 

PACS numbers: 72.80.Vp, 07.79.Cz, 72.10.Fk 


I. INTRODUCTION 

Graphene-based systems are promising candidates for 
the detection of the so-called bound states in the contin¬ 
uum (BICs)±£. BICs were first theoretically predicted 
by von Neumann and Wigner in 1929^ as quantum states 
with localized square-integrable wave functions, but hav¬ 
ing energies within the continuum of delocalized states. 
The electrons within BICs do not decay into the system 
continuum, thus these states should be invisible in trans¬ 
port experiments. 

The subject experienced revival after the work of Still- 
inger and Herrick in 19754 Since then, BICs were pre¬ 
dicted to appear in a variety of electronic, optical and 
photonic systems444 In these systems, effects of Fano 
interference^ were proposed as the underlying mechanism 
for the emergence of BICs and their possible experimen¬ 
tal observation. In particular, we recently proposed that 
BICs can be observed in the system of graphene with two 
collinear adatoms in the phantom atom configuration^. 

In this work, we show that the setup outlined in Fig{T| 
for suspended graphene can undergo a quantum phase 
transition (QPT) into the state with magnetic BICs if 
non-local graphene-adatom couplings are taken into ac¬ 
count. The phenomenon is a consequence of the partic¬ 
ular scaling of the local density of states (LDOS) T >0 on 
energy e in the vicinity of the Dirac point. The latter 
is proportional to the quantity known as pseudogap A, 
related to the intensity of the scattering near the Fermi 



Figure 1. (Color online) (a) Side view: two adatoms labeled 
by 1 (upper) and 2 (lower) placed collinear to a carbon atom 
beneath an STM tip in suspended graphene (inset) (b) Top 
view: the adatoms are coupled to the carbon atom beneath 
them and its nearest neighbors. The relative strength of these 
couplings define the Fano factor of interference qo playing the 
role of the natural control parameter of the system. It can be 
tuned by varying the slope of the Dirac cones in the graphene 
band-structure (inset). 


energy&4 Formation of the magnetic BICs becomes pos¬ 
sible only if A oc |e| 2 similar to the transition reported 
in Ref. [&] for a pair of quantum dots coupled to metallic 
leads. 

The QPT reported here is driven by a Fano factor 
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of interference qo which can be thus considered as the 
natural control parameter of the system. It can be 
tuned by changing the slope of the Dirac cones in the 
graphene band-structure (see FigtIJb)). The magnetic 
BICs appear within the region inside the critical bound¬ 
aries q c 1 < qo < q c 2 , where the dominant scaling law 
for the pseudogap is quadratic (A oc |e| 2 ). Outside this 
region, the mixed scaling A oc |e|,|e| 2 for qo < q c \ or 
cubic scaling A oc |e| 3 for q 0 > q c 2 , leads to the forma¬ 
tion of spin-degenerate BICs. The transition towards the 
magnetic BIC state is triggered due to the onset of the 
effective intrinsic ferromagnetic exchange field be¬ 

tween the adatoms mediated by the graphene monolayer. 


II. THE MODEL 


single carbon atom (local coupling regime). Otherwise, 
q 0 ^ 0 denotes the hybridization of the adatoms with the 
three second neighbors of carbons as depicted in Fig. [I] 
(non-local coupling). 

To analyze the transport properties of the geometry 
we consider and look for the existence of the BICs, we 
should focus on the local density of states of the host 
(LDOS) and those corresponding for the adatoms (DOS). 
The former defines the conductance of the device at zero 
temperature T = Oi: 

G ~ ^T tip LDOS, (3) 

with T t ip = 7rf 2 /?tip, Ptip as the STM tip density of states. 

To obtain the value of LDOS probed by the STM tip 
of Fig. [T] we should consider the tunneling Hamiltonian 


To give a theoretical description of the setup plot¬ 
ted in Fig[T] we use the model based on the Anderson 
Hamiltonia n 10 ' 11 : 


H 2 d = Y dk ^Fk)c\ ka c sk(J + ^2sjdd] a d ja 

sa ja 

+ ldJ2n djt n dA + J2 j dkV k(cl ka d ja +H.c.), 


(1) 


Tdtun — t c ^ ( t^tip.cr^ f(T "I" H.C., (4) 

<7 

where ^tip.o- and ’F 0 . are respectively fermionic operators 
for the edge site of the STM tip and 
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with vf being Fermi velocity. The graphene monolayer 
is described by operators c' skcr {c ska ) for creation (anni¬ 
hilation) of electrons in quantum states labeled by the 
wave number fc, spin cr and valley index s = 1,2. For the 
adatoms, dj a ( dj a ) creates (annihilates) an electron with 
spin a with energy Sjd, where j = 1,2 correspond to the 
upper and lower adatoms, respectively. The third term 
in Eq.([T]) accounts for the on-site Coulomb interaction 
U, with ridja = dj a dj a . Finally, the last term mixes the 
graphene and the levels £jd, wherein H.c. gives the Her- 
mitian conjugate of the first part. This mixing is charac¬ 
terized by the coupling V k = ^y/^-\Z\k\ v o(^-qo^g^), 
where A f is the number of conduction states, Ho denotes 
the unit cell area, and 


is the Fano factor of interference defined according to 
the results of Ref. [l2j . The parameter t stands for 
the coupling strength between carbon atoms, while vq 
and Vi represent the host-adatom hybridizations out¬ 
lined in Fig.[T| and D = 7eV denotes the band-edge for 
vf ~ c/300. The Fano factor q 0 can be tuned assisted by 
a variation of up, which enters into t = I^wp^ and mid. 
The experimental tuning of vf can be achieved, for in¬ 
stance, by means of modifying the carrier concentration 
in suspended graphene^ [inset of Fig. Oja)]. 

The situation q 0 = 0 corresponds to the scenario in 
which collinear adatoms are locally side-coupled to a 


is the field operator accounting for the quantum state 
of the graphene site placed right beneath the tip with 
hopping terms {td 1 and f c ), cf. Fig.[Q LDOS then can 
be computed as 

LDOS = ~~Y HGAe + )} = 2 V 0 + Y ALDOS, 7ff , 

7 r zz ' 

a ojl 

( 6 ) 

where C/ CT (e+) is the time Fourier transform of the Green’s 
function 

Sa = ~9 (r) Tr{p 2D [^(r), *+ (0)]+} (7) 

and 

= D5 (1 - ®I) )2 (8) 

is the graphene DOS, ALDOS^ stands for the part in¬ 
duced by the adatoms (see detailed derivation for it in 
the Appendix). 

It is worth mentioning that ALDOS^o- for j ^ l repre¬ 
sents electronic waves of a given spin a that travel forth 
and back between the upper and lower adatoms showed 
in Fig.[T](a), which for a specific energy e, become phase 
shifted by tt (Fano dip) with respect to the waves scat¬ 
tered by the adatoms, which are described by ALDOSjj CT . 
As discussed in Ref. [lj , such scattering process then pro¬ 
vides a mechanism of the emergence of BICs. This effect 
can be captured in the detailed derivation of LDOS ap¬ 
pearing in the Appendix. 
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According to the Appendix, the evaluation of 
ALDOSjjo- depends on the Green’s functions Gd ja d lrr 
(j = 1,2 and l = 1,2 ) for the adatoms. Additionally, 
to perceive the BICs emergence in our system, we should 
know the density of states DOST of these adatoms, which 
are determined as follows: 

DOS?- = (9) 

Thus both ALDOSjio- and DOST can be found by em¬ 
ploying the Hubbard I approach^ at T = 0, since 
the determined Hubbard bands match with those ob¬ 
tained via the Numerical Renormalization Group, in 
particular, for graphene with a single adatomi^. As 
a result, we can safely extrapolate the Hubbard I 
method to our graphene system. We start employing 
the equation-of-motion method to a single particle re¬ 
tarded Green’s function of an adatom in time domain 
Qd la i ja = ~i e ( T ) Tr{p 2 D [dz<x (t) , 4r (0)]+}, where d (r) 
is the Heaviside function, g 2 D is the density matrix 
of the system described by the Hamiltonian of Eq. m 
and [• • • ,•••]+ is the anticommutator between operators 
taken in the Heisenberg picture. Performing elementary 
algebra one obtains in the energy domain: 

( e "*~ — £ ld)Gdi a dj a = $lj + S 'y ' Gdj^djv 

i 

+ MGdi a n dlS ,d ja , (10) 


where e + = e + *0 + and 


s = £ / <*: 


v k v k 


■+ _ 


fivpk 


V 0 ft £ \2i 

-—e(l~q 0 -) In 


D 2 -e 2 


^q 0 (2- qo ^)-iA 


( 11 ) 


is the self-energy. Its imaginary part A is proportional 
to the scattering rate of the quasiparticles and is known 
as pseudogap. The latter is proportional to the local 
density of states of the host V 0 given by Eq.©, i.e., 
A = /TOo'Dcr'— . Thus 


Fourier transform, which leads to 


( £+ e Zd M)Gdi^n dl &,d ja — $lj < Tldit7 > 

+ dkVk @c ak „d* s d ls ,dj« 

S J 

_ n i n \ 

k <7 CT ’ dj <r d^-Cgkjjdifj^djfj 5 


(13) 


expressed in terms of new Green’s functions of the same 
order of Gd la n dlS ,d ja and the occupation number < n dld > 
determined by 


I r£ F —0 

< Ti dlS >=-/ Irn {Gd lB d ls )de. (14) 

7T J-D 


We highlight that for the quadratic pseudogap, the self- 
consistent evaluation of the Eq. (fl4l) reveals a range of 
magnetic solutions with < n dl -j- >^< n dl i > for the val¬ 
ues of go- lying in the range between two critical points 
q c i and q C 2 - Outside the magnetic region, i.e., for differ¬ 
ent scalings of the pseudogap (A oc |e|, |e| 2 for g 0 < g c i 
and A oc |e| 3 for g 0 > g C 2 ), Eq. (flT)) has non-magnetic 
solutions with < n d ,-[ >=< n dl i > only. This point will 
be addressed in detail in Sec.Ill of the paper (see in par¬ 
ticular Fig©. 

Furthermore, by employing the Hubbard I approxima¬ 
tion, we decouple the Green’s functions in the right-hand 
side of Eq. (1X51) as performed in Ref. [l| . This proce¬ 
dure enables us to solve the system of Green’s functions 

within Eq. (HU1), leading to Gd ia d, a = - 3 ■ where 

e-e Jd -^ - 


*? = (! + 


IA 'C.ti dj a 
e—ejd—U — Y. 


), and: 


£T=£ + AjAf-- 


(15) 


is the total self-energy, with j = 2,1 respectively for j = 

1,2 in order to identify distinct adatoms and Gd ia d- ja = 

>-^G d ~. d-. . , 

—-— Ja A* 7 are mixed Green s functions. 

£ Sj d ^ 


III. RESULTS AND DISCUSSION 


A = ^|£|d-^) 2 ’ (12) 

which depending on the value of the Fano parameter, the 
main contribution to the pseudogap can be linear, cubic 
or quadratic. As we will see in the discussion section, the 
latter situation is of particular interest, since magnetic 
BICs are formed in this case. 

In Eq. (fTOl) Gdi a n dl s,d ja is a t w0 particle Green’s func¬ 
tion composed by four fermionic operators, obtained 
from the time Fourier transform of Gdi a n d a ,d ja = 

_ i 6, ( T ) Tr {^ 2 D[dza (r) n dlS (r) ,d' jg . (0)] + }, with n dlS = 
djg-dia and spin a (opposite to a). Thus we first calcu¬ 
late the time derivative of Gd la n dlS ,d ja and then its time 


In the simulations we adopt T = 0 and the set of 
parameters^: £j d = £ d = —0.07D, which is feasible in 
suspended graphene (inset of Fig.QJa))^ and U = Vq = 
—2£ d - Additionally, to avoid that BICs decay into the 
continuum, we use t dl /t c = 0, otherwise it leads to ex¬ 
perimental detection of BICs by means of the so-called 
quasi-BICsi. 

In Fig© three distinct regions in the occupation num¬ 
bers of Eq. m for j — 1,2 appear identified by their cor¬ 
responding pseudogaps A [Eq. (fl2l) ]: the non-magnetic 
regions corresponding to small or big Fano factors appear 
to be divided by a magnetic central domain delimited by 
the critical values g c i and q c 2 . At critical values, abrupt 
jumps in the occupation numbers point out the existence 
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Figure 2. (Color online) Occupation numbers given by 
Eg. (1141) for spin-up and spin-down states of the adatoms as a 
function of go¬ 


of a QPT connected with the spin degree of freedom. 
Panel (a) of Fig. [3] presents the DOS corresponding to the 
regime qo = 0.8 < q c i where one can clearly see resolved 
and spin-degenerate peaks in Eq. © for the DOST. In 
Fig.©b) spin-polarized peaks emerge when the Fano fac¬ 
tor is placed within the boundaries q c \ < qo = 1.2 < q C 2 , 
while in panel (c) the case of qo = 2.0 corresponds to the 
limit of the phantom atom considered in detail in Ref. [T| 
for which spin degeneracy is recovered. 

To demonstrate that the system possesses BICs, we 
compare the density of states DOST for adatoms shown 
in Fig© with the host local density of states ALDOSy^ 
depicted at Fig© As one can see, both quantities reveal 
pronounced peaks (resonant states) placed at the same 
positions. Particularly in panels (a) and (b) of Fig. ©with 
qo = 0.8, we observe as aftermath of Ea. llA.13ll . degen¬ 
erate spin-up and down components for the Fano dips 
of ALDOSjjo- (l ^ j) interfering destructively with the 
peaks found in ALDOSyjo-. As this interference is com¬ 
pletely perfect, BICs emerge at the positions marked by 
vertical lines crossing panels (a), (b) and (c) of this figure. 
In panel (c) of the same figure, the total LDOS of Eq. © 
reveals absence of peaks at those places in which such a 
destructive interference occurs within panels (a) and (b). 
The aforementioned positions without peaks in Fig.©c) 
thereby give rise to BICs: the total LDOS that deter¬ 
mines the conductance does not catch the same peaks 
found in Fig. ©a) for the adatoms. Thus the aforemen¬ 
tioned invisibility of such resonant states points out that 
electrons with opposite spins stay equally trapped within 
these adatoms when qo < q c i and the pseudogap scales 
as A oc |e|, |£| 2 . 

Panels (a), (b) and (c) of Fig© correspond to the case 
<7ci < <Zo = 1-2 < q c 2 where magnetic solutions become 
possible, since the pseudogap is ruled by A oc |e| 2 . The 



e/D 


Figure 3. (Color online) (a) DOS for the case qo < q c i. Well 
resolved spin degenerate peaks are clearly visible, (b) DOS 
for the case q c i < qo < q c i with break of spin degeneracy, (c) 
DOS for the case qo > q c 2 when spin degeneracy is recovered. 


position of magnetic BICs is denoted by vertical dashed 
lines. Consequently, in the domain q c i < qo < q c 2 , the 
novelty due to a non-local coupling between graphene and 
collinear adatoms lies on the possibility of tuning the spin 
of the electrons trapped in the BICs of the adatoms. Such 
a feature yields an emerging based suspended graphene 
spintronics, in which a spin-filter of BICs rises as a fea¬ 
sible application. Outside the critical domain, just spin- 
degenerate BICs exist. 

Let us now present the physical arguments that eluci¬ 
date the emergence of the reported magnetic BICs, which 
is indeed triggered by a QPT. Similar QPT appearing due 
to the quadratic scaling of the pseudogap A oc |e| 2 and 
related breaking of the spin-degeneracy was discussed in 
Ref. 3; where a double dot system was explored. In re¬ 
gard of this dot setup, we highlight that the pseudogap 
A oc |e| 2 is only revealed to be present after performing 
a mapping of the original Hamiltonian into an effective 
model, in particular, under restricted constraints. On the 
other hand, we demonstrate that graphene emerges as 
the natural platform wherein the pseudogap of Ea. (fT2l) 
includes not only the regime |e| 2 , but also |e|,|£| 2 and 
|e| 3 , just due to the non-local adatom-graphene coupling. 
These regimes are accessible by means of the tuning of 
the Fano factor qo, which here is proposed to be prac¬ 
ticable by developing the Fermi velocity engineering^. 
Moreover, the non-local coupling assumption improves 
the emulation of the experimental reality, since the stan¬ 
dard case of local coupling regime, which is widely em¬ 
ployed in the literature, is indeed ideal and hides com¬ 
pletely the reported QPT. 

Fig ©illustrates how spin-resolved DOS of the adatoms 
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Figure 4. (Color online) Host local density of states corre¬ 
sponding to the cases of non-magnetic BICs (panels (a), (b), 
and (c)). BICs appear when a peak in Ej ALDOSj J(T is fully 
compensated by a Fano dip in Ej^j ALDOSjjo-. Positions of 
BICs are marked by vertical dashed lines. Panels (a) and (b) 
correspond to spin resolved ALDOS. Lower panel (c) corre¬ 
sponds to total LDOS defining the conductance. 


depend on the Fano parameter. The variation of go in the 
wide range below the critical value g c i ~ 1.1766 shifts the 
position of the peaks corresponding to opposite spin com¬ 
ponents equally, as it is shown in the upper panel. How¬ 
ever, above the critical value the spin splitting abruptly 
appears as it is shown at the lower panel of the figure, 
which clearly indicates that the system undergoes a QPT. 
The abrupt appearance of the spin splitting is intimately 
connected with the step-like behavior observed in the oc¬ 
cupation numbers shown in Fig. [2] Note that the increase 
of the Fano factor above the second critical value q c 2 leads 
to the recovering of the spin-degeneracy as the regime of 
the phantom atom with cubic scaling of the pseudogap 
A ex |e| 3 is achieved. 

Within the critical boundaries q c ± < qo < q C 2 , the 
quantity R,e(£A — £) = J exch from Eq. (fl5l) plays the 
role of a Zeeman-like splitting of the levels Ed in the 
adatoms. This splitting arises from an intrinsic exchange 
field J rexch between the adatoms intermediated by the 
graphene monolayer. Its value is ruled by the system nat¬ 
ural control parameter, namely the Fano factor go, which 
drives the graphene system towards a QPT. As the upper 
and lower adatoms magnetize equally, cf. Fig.[2j the cou¬ 
pling between them is revealed as ferromagnetic. Note 
that the dependence of the effective held on the Fano 
parameter is non-monotonous: it drops abruptly when 
go = Qc 2 ~ 1.3582. 


Figure 5. (Color online) Host local density of states corre¬ 
sponding to the cases of magnetic BICs (panels (a), (b), and 

(c))- 
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Figure 6. (Color online) (a) Spin-degenerate crossover from 
merged peaks of Eq.© for the adatoms DOS towards to re¬ 
solved peaks, (b) QPT due to an abrupt spin-splitting of the 
peaks. 


IV. CONCLUSIONS 

In summary, we have proposed a setup based on 
graphene-adatom system in which magnetic BICs are 
triggered by a quantum phase transition in the region 
of the quadratic scaling of the pseudogap with energy, 
A x |e| 2 . The control parameter which drives this tran¬ 
sition is a Fano factor of interference tunable by changing 
the slope of the Dirac cones in graphene band-structure. 
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Appendix: LDOS derivation 

To obtain the analytical expressions of the LDOS 
given by Eq. © appearing in the conductance of Eq. © , 
we begin by applying the equation-of-motion approach 
to g a = -j-9 (r) Tr{£» 2 D[^ CT (T), 4'J(0)]+}, with Eq.© 
rewritten as 

+ (ttVqVo) E Cj dj Q -, (A . 1) 

3 

expressed in terms of Cj = {'kVqVq)~ 1 (t^ / t c )5ji . Substi¬ 
tuting Eq. (|A.1I) in one finds 

X VWli 1 - 9o r ^-)dqQc eka c- sq<r + (ttDo^o) 

js 

X {Gd ja c ska + Gc ska d ja ) + ( n’Do'Vo ) 'y , B jClGd ilT di^, 

jl 

(A.2) 


with the new Green’s functions Q r t r - , 
Sd ja c Bka , Qc ska d ja and Qd ja d to be deter¬ 
mined. To this end, we first consider Q r , r - = 

-jrd (r) Tr{p 2 D[c s fca (r), c\ q<7 (0)]+}, whose time deriva¬ 
tive d T = gives 

9rGc Bkq c Bqq = ~^S (T)Tr{g 2D [c ska (r), c \ qa (0)] + } 

— Jj^ VF ^^ CskcrCM i'’ _ Ji ^ ' J VkGd jq Csq<,i 

3 

(A.3) 

where we have used 


ihd T c s ka (r) = [c sfco .,'H 2 d] = ( hv F k)c sk(J (r) 

+ E Vkd i°- ( T ) ■ 

3 


In the energy domain after performing the time Fourier 
transform, we solve Eq. (1A.3I) for Gc Bka c Bq a and obtain 


- _ 6 (k — q) 5 S § V fc p 

yc ek „c Sq< , £ + _ hvpk + 2^ £ + _ h,VFk ddi '’ Csqa ' 

(A.5) 

Notice that we also need to calculate the mixed Green’s 
function Gd ja c aqa ■ We then define the advanced Green’s 
function Td jqCaqq = (~t) Tr{^ 2D [d] CT (0) , c sq(T (t)]+}, 

whose equation-of-motion reads, 

dTFd ja cs q a =~\ 8 ( r ) Tr (92D \d) a (0), c sqa (r)]+} 

^ ~jl^ VF ^ d ~ d j' TCa '> Cr ^ ^ y '^q-^djcdlg i (A.6) 

i 

where we have used once again Eq. (IA.4I) . interchanging 
k «->• q. The Fourier transform of Eq. (IA.6I) leads to 

e J~djc,c sq a — (hv F q)fd jqCsq „ + y ' Vq-Fdivdiai (A.7) 

i 

with e**=e- *0+. Applying the property Gd ja c sq „ = 
{Fd jqCaqq y on Eq. (1A.7D , we show that 

E+Gdj q c Bq<7 = (fcVFq)Gdj q c sq<7 + y ', VqGd icr di a i (A.8) 


Qd jaCaqa E £ + _ hvpq 


^ djcr diu 


(A.9) 


and analogously, 


^CsqcT djc 


= E 


v 0 


£ + — hv F q 


Gdi^dja ■ 


(A.10) 


Now we substitute Eq. (1A.9I) into Eq. (1A.5I) and the 
latter, together with Eqs. (1A. 101) and (TTTT) for the self¬ 
energy splitted as 


E = E / dk £ + V - V n VF k = ^oiAj - iBj), (A.11) 

into Eq. (IA.2I) in the energy domain, which results in 



+ (ttV 0 v 0 ) 2 - iB 3)Gd ja d,A-Ai - iBi) 

31 


+ (7rP 0 t , 0 ) 2 y ^Cj(Ai - iBi){Gd ja d iq + Gdivdjq) 

ji 

+ ( ttT> 0 vo ) 2 ^ CjCiGd^div- 
31 


(A.4) 


(A.12) 







































Thus after some algebra via the evaluation of 
— i I m (^cr): we determine Eq.© as the LDOS probed 
by the STM tip, with 
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ALDOS, 7ff = -(nv 2 0 Vl)lm[(Ai - iBi)Q dladja (Aj - iB 0 )\, 

(A-13) 


with Aj 
1 


= 1/2 (tdi/t c ) and#,- = 

ImE. 
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